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A boojum is a topological defect that can form only on the surface of an ordered medium 
such as superfluid "^He and liquid crystals. We study theoretically boojums appearing between 
two phases with different vortex structures in two-component BECs where the intracomponent 
interaction is repulsive in one phase and attractive in the other. The detailed structure of the 
boojums is revealed by investigating its density distribution, effective superflow vorticity and 
pseudospin texture. 
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The study of topological objects is important because 
such objects appear not only in condensed matter physics 
but also in cosmology and high-energy physics.^ Atomic- 
gas Bose-Einstein condensates (BECs) are ideal for ex- 
amining topological defects and textures in a quantum 
condensed system.^ A major advantage of this system 
is that it can be accurately described using the mean- 
field theory. Moreover, it is a good system to study ex- 
perimentally because the atomic interaction is tunable 
through the Feshbach resonance and optical techniques 
allow one to control and visualize the condensates di- 
rectly. Although topological defects and textures have 
been studied thoroughly in anisotropic superfluid •^He,^-'' 
interest in BECs with multicomponent order parameters 
has been increasing. For example, Kasamatsu et al? nu- 
merically simulated vortex structures in two-component 
BECs that were similar to those in superfluid ■^He-A. 
Hence, we can learn about topological defects and tex- 
tures in multicomponent BECs as well as in superfluid 
^He. This sort of work can lead to remarkable develop- 
ments toward the elucidation of topological objects in 
nature. In this study, we theoretically examine boojum 
structures in two-component BECs and show how they 
are analogous to boojums on the interface between A 
and B phases of superfluid '^He. This study pioneers the 
research field of interfacial topological defects of multi- 
component BECs. 

A boojum is a point defect that can exist only on the 
surface of an ordered medium. The name came from Mer- 
min who analyzed these defects in superfluid "^He.^ Boo- 
jums at the interface separating A and B phases of super- 
fluid ■^He have been recently studied by Blaauwgeers et 
al.^ They made the A-B phase boundary in a cylindrical 
container by controlling the gradient of the magnetic field 
along the axis of the cylinder. When the container was 
rotated, each phase had a distinct structure of quantized 
vortices. In particular, the vortex core diameter of the A 
phase is 10"^-fold larger than that of the B phase. When 
the angular velocity Q is smaller than a critical value flc, 
the A phase rotates by making doubly quantized vor- 
tices with "soft cores", whereas the B phase does not 
rotate. Also, the A phase vortices bend near the phase 
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Fig. 1. Boojums in rotating two-component BECs. (a) The con- 
stant surface density + |V'2p is shown as the blue surface. 
The vortices of component 1 are red and those of component 
2 are light green, (b) Vortices of component 1 (left) and com- 
ponent 2 (center), and superposition of two vortices(right). The 
gray planes are the z = planes, (c-e) Density distributions 
(left) and |V'2p (right) in various cross sections perpendic- 
ular to z. (c) Section at 2 = 10.8. (d) Section at z = —0.4. (e) 
Section at z = —11.2. The box dimensions along the x, y, and z 
axes are -12.8 to -1-12.8, -12.8 to +12.8, and -25.6 to -1-25.6, 
respectively (256 X 256 X 512 in discrctized space). 



boundary toward the container wall.^ In contrast, when 
n exceeds flc, the A phase vortices penetrate into the B 
phase, which can also rotate by forming single-quantized 
vortices with " hard cores" . The vorticity distribution is 
continuous in the A phase, but singular at the cores of 
the vortices in the B phase. The boojums appear at the 
A-B phase boundary where the two kinds of vortices con- 
nect. However, we have little quantitative understanding 
of the structure of the order parameters around boojums. 
This lack of knowledge is partly a result of the difficulty 
in analyzing a system with such a large difference in co- 
herence length between the two phases. 

Because of the similarity of the order parameters be- 
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tween superfluid ^He and a liquid crystal, boojums have 
also been studied in liquid crystals.^' ^ Boojums are also 
reported in molecular patterns of fatty acids at an air- 
water interface.^" Then, atomic BEC, which is a well- 
controlled quantum system, should be a ground giving 
boojums. The research in this direction may push out 
multicomponcnt BECs as a testing ground of the physics 
of pattern formation in which various topological defects 
play important roles. ^ 

For this study, we numerically created two phases of 
different vortex structures in two-component BECs elon- 
gated along the rotation axis n = Qz, as shown in Fig. 
1(a). The intercomponent interaction is tuned to be re- 
pulsive for 2: > and attractive for z < 0. Thus, for 
repulsive interactions, the vortices of one component al- 
ternate with those of the other component, whereas the 
vortices of two components have common cores for at- 
tractive interactions, -^^"^^ as shown in Figs. 1(c) and 1(e). 
Then, just below the z = plane, bifurcation points of 
vortices form (Fig. 1(b)), which create boojums. The re- 
gions above and below the boojums are called the a- 
and /3-phases, respectively. By introducing an effective 
velocity and the pseudospin representation, we can more 
easily examine the structure of boojums and reveal their 
relationship to textures in superfluid ^He. 

We obtained stable states of two-component BECs in a 
frame rotating with an angular velocity ft = Qz by solv- 
ing the time-independent, coupled Gross-Pitaevskii (GP) 
equations. For simplicity, we assume that m = mi = m2, 
N = Ni = N2, and u!± = uji± — uj2±, where rrii, Ni and 
LOi± are the particle mass, the number of particles, and 
the radial trapping frequency of component i {i = 1,2), 
respectively. It is convenient to measure the length, time, 
and energy scale in units of b± =\fhJrm}Z^ wj, and /k<;_L, 
respectively. In these units, the coupled GP equations are 



(i = l,2), (1) 

where pn is the chemical potential of component i, Lz = 
—i{xdy — ydx) is the angular momentum operator, and 
the wave functions are normalized as j dr\ipi\'^ = 1. 
The harmonic potential is V = ^ix^ + + >?z'^) with 
anisotropic parameter A = 0.3. The intracomponent 
and intercomponent coupling constants respectively are 
C^ii = U22 = 4:TTNau/bj_ and U12 = U21 = AnNau/hi, 
with the corresponding s-wave scattering lengths an and 
ai2 being dependent on z. 

In this system, it is essential for making boojums that 
BECs have the two regions with U\2 positive and nega- 
tive. Then we assume U\\ = U22 = 1000 and U12 = 750 
for z>0, Un = U22 = 2500 and U12 = -750 for z < 0, 
and 12 = 0.8. For making stable boojum structures, it 
is neither necessary to change the interaction parame- 
ters so sharply nor to make Un and U22 spatially de- 
pendent. Such a parameter setting is just for visualizing 
the structures of boojums clearly. The density distribu- 
tions of each component are shown in Figs. l(e)-l(j). 
The choices of D, and U12 were based on the phase dia- 
gram by Kasamatsu et al}^ such that the a-phase vor- 



tices made the alternating square lattice that is shown in 
Fig. 1(c). In the /?- phase, the two components have the 
same density distribution; that is, ni(r) = 7x2 (r) where 
ni{r) = \ipi{r)\^. In this case, the two components are re- 
duced to one component, thus making a triangular vortex 
lattice (see Fig. 1(e)) with a solid-body rotation except 
near the vortex cores with singularities. Thus, the vor- 
tices of each component are bifurcated on the interface 
between the two phases. The structure of one pair of two 
vortices is shown in Fig. 1(b). Here, a boojum exists at 
the point where the pair of two vortices is bifurcated. 

The structures of the boojums of this system were re- 
vealed by investigating the distribution of the vorticity 
of the effective superflow velocity VeS- This velocity is 
defined as^^ 

VeS = : , (2) 

ni -I- 712 

where = riiWOi and 9i is the phase of ipi. From Vcff, 
we obtain the effective vorticity as (jJcS = V x Vcs- 
The resulting three-dimensional structure of |a;off| near 
the boojums is shown in Fig. 2(a), which may be called 
a boojum lattice. Figures 2(b)-2(d) show the distribu- 
tions of the z component {(jJes)z and the unit vector 
Veff/l'WeffI in three cross sections perpendicular to z. In 
the /3-phase, the vorticity distribution is similar to that 
of a one-component BEC with singularities at the vor- 
tex cores (Fig. 2(d)). We call such a singularity of o^eff 
"singular vortex". In contrast, in the a-phase, the vor- 
ticity distribution is continuous, making a square lattice 
of fuzzy round blobs (Fig. 2(b)). Each blob has the vor- 
tex core of either component (red and blue points in Fig. 
2(b)). We call such a blob "continuous vortex". A boo- 
jum appears at the point on the interface where a singu- 
lar vortex splits into two continuous vortices. 

The bifurcation of vorticity through the boojum was 
found by investigating the circulation. The circulation 
per singular vortex in the /3-phase is 27r, just like that in 
a typical one-component BEC. The circulation around a 
continuous vortex was estimated by first assuming that 
the square lattice is periodic. Then, each Wigner-Seitz 
cell has a continuous vortex, and ni = n2 on the bound- 
ary of the cell.^^ The circulation per continuous vortex 
is 

/• ,„ / niV^i+n2V^2 
/ Weff ■ db = j • dl 

J cell ^boundary f^l + ^^2 

= J / (V^l + V^2) ■dl=7T. (3) 

^ J boundary 

Thus, the whole circulation per singular vortex in the /3- 
phase is transmitted to a pair of continuous vortices in 
the a-phase. This process of vorticity bifurcation makes 
the velocity field in the a-phase more similar to a solid- 
body rotation than that in the /?-phase because the num- 
ber of vortices of each component in the a-phase is al- 
most equal to that in the /3-phase. This is observed by 
comparing the velocity fields in the a-phase without sin- 
gularities (Fig. 2(b)) and in the /3-phase with singulari- 
ties (Fig. 2(d)). 

This consideration makes us think of the interest- 
ing structure of boojums of an n-component BEC; a 
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Fig. 2. Boojums represented by vorticity of effective superflow. 
(a) Three-dimensional distribution of |ttiofi | with several boo- 
jums near interface, (b-d) Distributions of (tti^ff )z and vectorial 
plot of unit vector I'off / 1 1'cff I in cross sections perpendicular to z 
at (b) 2 = 10.8, (c) 2 = -0.4, and (d) z = -11.2. Red shows the 
region with high vorticity. Plot (b) shows the lattice of continu- 
ous vortices and four Wigner-Seitz cells in the o-phase. Red and 
light blue points in the cells show the positions of the vortices 
of components 1 and 2, respectively. Plot (c) shows the pairing 
of continuous vortices near the interface. Plot (d) shows that 
singular vortices make a triangular lattice below the boojums 
(/3-phase). 



singular vortex may split at a boojum (n-boojum) to 
n continuous vortices with the same circulation 27r/n. 
When n — > cx), vorticity becomes uniform {v^s ^ 'U solid, 
(jjcs — > 2J1) above the lattice of oo-boojums, as shown in 
Fig. 3(a). Such a topological structure of vorticity with 
oo-boojums should generally appear on a boundary be- 
tween one medium with singular vorticity and the other 
with a uniform vorticity. One example is the interface 
between a type-II superconductor and a normal media 
under a magnetic field perpendicular to the interface, 
where oo-boojums of magnetic flux should appear. 

It is useful to understand the topological structure of 
vorticity of a boojum to consider its analogy to a Dirac 
monopole. When we consider the vorticity Weff as the 
magnetic flux density B, a Dirac monopole and a string 
respectively correspond to a boojum and a singular vor- 
tex.^ However, the radiation of uictf from a boojum is dif- 
ferent from that of B from a Dirac monopole, as shown 
in Figs. 3(b) and 3(c). 

The pseudospin representation^'* provides insights for 
understanding the similarity between two-component 
BEG and superfluid ^He. We consider a spin-1/2 BEG 
for the following spinor ^ = [ipi,tp2]'^ = i'rX ^ tti^ 
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Fig. 3. (a) Lattice of oo-boojums. (b) Vorticity distribution 
around boojum. (c) Magnetic flux distribution around Dirac 
monopole. 
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X = 



cos(i?/2)e-*'^/2 
sin(t?/2)e*'^/2 



(4) 
(5) 



ni+n2, 9t = ^1 + ^2, cosi? = (ni-n2)/(n-i + 
= 02 ~ Oi. The parameters and ip refer to 



where tit = 
712), and ip 

the polar and azimuthal angles of the local pseudospin s 
which is defined as s = X^^X = cos?? 2: -I- sin'!9(cosiy9a; -|- 
siniysy), where cr is the Pauli matrix and = 1. From 
these relations, the effective velocity Vcs of eq. (2) and 
its vorticity Woff are 

1 ^v*^ 



■Wcff 



2i 



icos??V<y9, (6) 



-icosi?V xVlp+ ^(Vi?) X {sm-dVif). (7) 

The first and second terms in eq. (7) vanish except at the 
vortex cores. The scalar product [(V??) x (sim? V(,c)] • dS 
is equal to the infinitesimal solid angle dfl covered by 
the s orientations within the infinitesimal plane dS. As 
a result, the third term can be rewritten as 

i(Vi?) X (sindVp) = ^e,JkS^VsJ x Vs^ = q{r), (8) 

which is connected to the Mermin-Ho relation for the s 
texture.*^ Figure 4 shows the pseudospin texture and the 
values of ip in cross sections perpendicular to z. In the 
a-phase (Fig. 4(a)), the first and second terms cancel, 
leading to a; off — q{r). Thus, the circulation of cq.(3) re- 
quires that all points on the surface of a unit hemisphere, 
drawn with s as the radius vector, are covered by the s 
orientations within a Wigner-Seitz cell in Fig. 2(b), as 
observed in Fig. 4(d). The pseudospin texture in the a- 
phase is similar to the texture of the order parameter I in 
the structure of a locked vortex (LVl) in a rotating su- 
perfluid '^He-A.^^ When we approach the interface, spin 
domains are formed with a domain wall that connects 
two pairs of up-spin and down-spin (pairs of the origi- 
nal vortices), as shown in Figs. 4(b) and 4(e). A boojum 
appears at the point where the up-spin and down-spin 
annihilate. This annihilation removes domain walls, thus 
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Fig. 4. Textures of pseudospin s and distribution of ip (the az- 
imuthal angles of s) in cross sections perpendicular to z. (a) 
and (d) Cross sections at z=10.8. (b) and (e) Cross sections at 
2 = —0.4. (c) Cross section at z = —11.2. Red and blue points 
show the positions of down-spin and up-spin at the vortex cores 
of components 1 and 2, respectively, (d) and (e) are magnified 
images of (a) (rotated by 45°) and (b), respectively. The domain 
walls are shown as pink dotted lines in (b). The four squares in 
(d) are the four Wigncr-Seitz cells in Fig. 2(b). 

making the domains larger. In the /3-phase, the spin do- 
main waUs vanish, which means that s becomes uniform 
and Ues in the x-y plane; x — const. Thus, the /3-phase is 
described only by ipx, which behaves similar to the order 
parameter of a one-component BEC such as superfluid 
3He-B. 

In conclusion, we have numerically simulated boojums 
in two-component atomic BECs by making the intercom- 
ponent and intracomponent coupling constants spatially 
dependent. The structure of the boojums is analyzed by 
introducing effective velocity and pseudospin representa- 
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tion. In general, boojums should appear in rotating two- 
component BECs with the intercomponent interaction 
repulsive in one phase and attractive in the other. It is 
experimentally possible to control the interaction param- 
eters by the Feshbach resonance tuned with a spatially 
dependent external field. Hence, we expect that boojums 
could be observed in two-component BECs. Studies on 
the dynamics of this system and BEC systems of more 
than two components are now in progress and will be 
reported elsewhere. 
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